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Abstract 

I present analytical expressions for the collinear and soft gluon corrections to Higgs 
production via the process bb ^ H as well as gg ^ H through next-to-next-to-next-to- 
leading order (NNNLO). The soft corrections are complete while the collinear corrections 
include leading and some subleading logarithms. Numerical results at the Tevatron and 
the LHC are presented, primarily for bb H. It is shown that the collinear terms greatly 
improve the soft and virtual approximation at next-to-leading order (NLO) and next- 
to-next-to-leading order (NNLO), especially when subleading terms are included. The 
NNNLO collinear and soft corrections provide significant enhancements to the total cross 
section. I also provide expressions for the collinear and soft corrections through NNNLO 
for the related Drell-Yan process. 



1 Introduction 



The search for the Higgs boson is one of the most important goals at the Tevatron and the 
LHC colhders p]. The main Standard Model production channel at these colliders is gg — > H. 
However, the channel bb H can also be important in the Minimal Supersymmetric Standard 
Model (MSSM) and other theories beyond the Standard Model. In fact, in the MSSM bb ^ H 
dominates gg —>■ H at high tan f5, where tan f3 is the ratio of the vacuum expectation values 
for the two Higgs doublets. Complete results through NNLO are known for gg ^ H [31 IH [3], 
bb ^ H ^ and the related Drell-Yan processes qq —>■ V [7|. 

Hard-scattering cross sections can receive large corrections from soft-gluon contributions 
near threshold [H [U [10]. Such corrections have been studied for a wide class of cross sections 
(for a review see [11]) where the soft corrections dominate the cross section. For the process 
gg — >■ H, however, it was shown in [12] that the soft (or even the soft and virtual) corrections are 
inadequate to serve as a good approximation of the complete cross section at NLO. Collinear 
corrections have to be added to obtain a reasonable approximation. This was also shown at 
NNLO in [5]. We verify these results for gg —>■ H, and we investigate the soft and collinear 
approximations for the process bb — > H. We find that at both NLO and NNLO the soft, or 
soft and virtual, approximation is inadequate and that collinear logarithms must be added to 
improve the approximation. 

The processes qq ^ V, bb H, and gg H, present a unique opportunity to study com- 
plete soft-gluon corrections even at NNNLO. Although partial NNNLO results have appeared 
for other processes, such as top quark production [131 E] ) the complications involved in such 
processes are far greater. The processes qq V, bb ^ H, and gg H, are much more 
amenable and easy to study (relative to, for example, heavy quark or jet production) for two 
distinct reasons: 

(i) The color structure of these processes is trivial. At lowest order the processes are 
electroweak and the final state does not involve any colored particles (quarks or gluons) so the 
color vertex is trivial. This is to be contrasted with processes such as heavy quark [L^ or jet 
[15] production where the complicated color structure requires matrices in the study of the soft 
corrections. 

(ii) The kinematics of these processes in the calculation of total cross sections are trivial. 
Again this is to be contrasted with heavy quark or jet production where the kinematics is more 
complicated and where the structure and form of the soft-gluon terms depends crucially on the 
kinematics. 

It should thus come as no surprise that Drell-Yan and Higgs production are more well 
understood and known than other hard scattering processes (save deep inelastic scattering). 

A unified approach to calculating the soft-gluon corrections for arbitrary hard-scattering 
processes through NNNLO is available in [13] . We note that related processes involving the 
Higgs to which the formalism has been applied are charged Higgs production via bg —>■ tH~ 
[131 [TB] and the process bg —>■ bH [T7]. The NNNLO master formula in Ref. is quite long 
because it is structured so it can address all the possible complications in color structure and 
kinematics of a general process. Because of the two reasons mentioned above the application 
of the formula to qq V, bb ^ H, and gg H, results in much simpler expressions. 

For the process bb ^ H the complete NNLO corrections were presented in [6j; however 
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the analytic results presented there did not include the scale-dependent terms, and the color 
structure of the corrections (which is useful to know from a theoretical standpoint and for 
resummation) was not made explicit. In this paper I present the complete soft-gluon cor- 
rections and the coUinear logarithms beyond leading accuracy at NLO, NNLO, and NNNLO 
including the full color factors and the scale-dependent terms. Analytic expressions for the 
scale-independent terms in the soft-gluon corrections through NNNLO were also presented in 
Ref. [18] including the full color structure. The results in this paper are in complete agreement 
with [18]. 

Studies have been performed on the soft-gluon corrections to the process gg ^ H as well 
as the Drell-Yan process at NNNLO in [181 IlH] • However the analytical results for the scale- 
dependent terms were not presented. In this paper we derive the full NNNLO soft corrections 
for gg H and qq — > V , including explicitly terms with the factorization and renormalization 
scales, thus verifying and extending the results presented in [TH [19] . Furthermore, we calculate 
coUinear logarithms through NNNLO. 

The goals of this paper can be summarized as follows: 

(1) For the process hb —>■ H: To extend earlier results for the NNLO soft and coUinear 
corrections; to investigate the contribution of soft, virtual, and coUinear corrections to the total 
cross section and determine how well they approximate the complete corrections at NLO and 
NNLO; to derive at NNNLO the complete soft-gluon corrections and the coUinear corrections 
beyond leading logarithm accuracy; and to provide for the first time a detailed study through 
NNNLO of the numerical impact of these corrections at the Tevatron and the LHC 

(2) To provide the complete analytical results of the soft-gluon corrections, including scale- 
dependent terms, and the coUinear corrections beyond leading accuracy through NNNLO for 
the process gg ^ H and for Drell-Yan production. Since these processes have been studied 
extensively before we do not study numerical applications of the results (except in a limited 
way for gg ^ H for comparison with hh —>■ H). 

(3) To study the effects of subleading terms in the soft-gluon expansion and in coUinear 
corrections and make connections with earlier results for various processes. 

The paper is organized as follows. In Section 2 we provide a short overview of the formalism 
for the resummation of soft and coUinear logarithms in Higgs and Drell-Yan production. In 
Section 3 we derive explicit analytical expressions for the soft and coUinear corrections for 
the process bb H. Detailed numerical results are provided in Section 4 for bb H a.t 
the Tevatron and the LHC, together with a few results for gg — > H. The conclusions are 
in Section 5. Appendix A collects detailed expressions for the numerous quantities needed 
in the calculations. Appendices B and C contain explicit analytical expressions for the soft 
and coUinear corrections for the Drell-Yan process and for Higgs production via gg —>■ H, 
respectively. 

2 Resummation of soft and coUinear logarithms 

It is well known that soft and coUinear logarithms in hard-scattering cross sections exponen- 
tiate [H [H [To]. In this section we review the resummation formalism for these corrections. 
The resummation of soft and coUinear logarithms is carried out in moment space. We define 
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moments of the partonic cross section by 6"(A^) = / dz z^~^d{z), with z = M'^/s. Here M is a 
hard scale which we will take to be the Higgs mass mn for Higgs production and the invariant 
mass a/Q^ of the produced dilepton pair for Drell-Yan production, while s = {pi +^2)^ with pi 
the momenta of the incoming partons. At partonic threshold, z = 1. N is the moment variable 
and the logarithms of exponentiate. The resummed partonic cross section in moment space 
can then be written as 



X exp 



exp 



-2 / dz— I I —A, {asiXM^)) + Vi - zfW 



4 / — liji 

J pup /i 



\- z 
exp 
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(2.1) 



Factorization-derived resummation studies are based on the above formula or its equivalents. 
In particular, resummation studies for simple-color-structure processes such as Drell-Yan and 
Higgs production are based on [HI [9] while processes with more complicated color structure 
require the additional ingredients in [IHl [HI [131 [IS] • 

The expressions for the quantities in the exponents up to three loops are given in Appendix 
A together with the necessary integrals for NNNLO calculations. The first exponent includes 
corrections that are universal in hard scattering cross sections and only depends on the identity 
of the incoming partons (quarks or gluons) [H [9]. The second and third exponents control 
the factorization scale, fip, and the renormalization scale, fiR, dependence of the cross section, 
respectively. H is the hard-scattering function for the scattering of partons, while S is the 
soft function describing noncollinear soft gluon emission [lOJ. Also A^ = Ne'^'^ where 7^; is 
the Euler constant. The evolution of the soft function follows from its renormalization group 
properties and is given in terms of the soft anomalous dimension Ts which can be explicitly 
derived through the calculation of eikonal vertex corrections. 

To include collinear singularities we essentially replace —2{z^~^ — — z) in the first 

exponent of Eq. (12.11) by 2z^~^. These corrections are of the form In"^ N/N in moment 
space (see Appendix A). The leading collinear corrections come from the one- loop term in Ai 
[T2I [20] . However at NNLO and beyond one does not derive all the subleading terms from this 
exponential. We will show, nevertheless, that both analytically and numerically we get the 
dominant terms and the approximation is excellent. 

The exponentials in the resummed cross section can be expanded to any fixed order in the 
strong coupling and then inverted to momentum space to provide explicit results for the 
higher-order corrections. A fixed-order expansion avoids the problems with infrared singularities 
in the exponents and thus no prescription is needed to deal with these in our approach (see 
discussion in Ref. [211 )• The n-th order corrections can then be written as 



2n-l 



5(1 -Z)+ J2 
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Here V^"^ are the virtual contributions, the soft contributions are of the form of plus distribu- 
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tions with coefficients 5*^"^ and C]l'' are the coefficients of the colhnear logarithmic contribu- 
tions. 

To calculate the hadronic cross section we convolute the partonic cross section with parton 
distribution functions 0. The n-th order corrections to the hadronic cross section can then be 
written as 

= T I rlT, I rl^^ I rio' S(t - rr.,^^^ M^A rh(^„\ o'\ (2.3) 



dxi / dx2 I dz 6{t — X1X2Z) (f){xi) (f){x2) d'^^\z) 
t/xi Jo 



Here r = Q'^/S and z = Q'^/s for Drell-Yan production, with S the center-of-mass energy 
squared of the incoming hadrons, while r = rn^j/S and z = m\/s for Higgs production. 
Substituting the expression of Eq. (12.21) in Eq. (12. 3p we have 
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which, after a few manipulations, gives 
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3 Soft and collinear corrections for hh ^ H through NNNLO 

In this section we calculate the complete soft-gluon corrections and a large class of collinear 
corrections beyond leading accuracy to the cross section for the process bb ^ H through 
NNNLO. In the calculation we also need the virtual corrections through NNLO [6| IT8]. 
We use the notation 

"inVl - z] 



Diiz) 



to denote the plus distributions in the soft corrections. 

The complete NLO soft and virtual corrections are then given by 
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where the color factors Cp-, Ca and the Q constants are defined in Appendix A, and the Born 
term, F^^^ is given by 

(3.3) 



For the Standard Model process Xbi^ji) — ^f^bilJ^'u) / v , where v = 246 GeV is the Higgs 
boson vacuum expectation value and rub^fiji) is the MS bottom quark running mass at scale 
fiji. In the MSSM the value of A;, depends on which supersymmetric Higgs boson is produced. 
In the numerical results below we will mostly show ratios of cross sections so the results are 
equally valid for the Standard Model and the MSSM. However, when values for cross sections 
are presented in this paper they are always for the Standard Model process. 
The leading and next-to-leading collinear corrections at NLO are 
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The complete NNLO soft and virtual corrections are given by 
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where is the number of light quark flavors. This is in agreement with Refs. [HI ITS] . 
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It is interesting to note here that the scale-independent contribution in the Dq terms is solely 

^(2) 



due to terms and the two-loop function Ggg defined in Appendix A. It was shown in Ref. [22] 

(2) 

that for top quark production the inclusion of Q and two-loop Ggg terms in the Dq coefficient 
for that process provided the majority of the Dq corrections, even though the complete two-loop 
Dq terms are not known for that process. This was demonstrated by comparing the corrections 
in two different kinematics formulations for ti production and showing that the results agree 
only with the inclusion of these Q and two-loop terms. 

The leading and some subleading coUinear corrections at NNLO are 
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Several remarks are in order here. The leading coUinear (LC) logarithms at NNLO (i.e. 
ln^(l — z)) are complete. The next-to-leading coUinear (NLC) logarithms (i.e ln^(l — z)) are 
not complete. However, numerically they are an excellent approximation to the complete NLC 
terms. Also analytically the CfCa-, ujCf and the ln(/i|n/m^) terms are exact. The next-to- 
next-to-leading coUinear (NNLC) logarithms (i.e ln(l — z)) are also not complete. However, 
again numerically they are an excellent approximation to the complete NNLC terms. Also an- 
alytically the CfCaC2i ln^(/^F/"^li-)) and ln(yu|j/m|^) terms are exact. Finally, in the con- 
stant terms we only show the scale-dependent terms; the ln^(/i|,/m^), ln(/i|./m^) ln(/i^/m^) 
and ln(/i|./m|^) terms are exact, while the \n{ii\/m'\j) terms are almost exact. More details 
on the numerical approximation are given in Section 4, and on the analytical approximation in 
the context of the related Drell-Yan process, where explicit results involving the color factors 
have been published [7J, in Appendix B. 

The complete NNNLO soft-gluon corrections are given by 



a 



{3)S 
bb^H 



+ 



8ClV,(z) 



-8CUi + GC2) + C'pCa 




+ 



110 
121 



+ —Clnf-20Cl In 



9 



m 



H , 



VAz) 



40 



44 



27 



27 



IQCl In^ 



H , 



-UCl+^-^ClCA-^-^Cln 



In 



m 



H , 



7 



+ ( 12C| 



44 



+ 



169 



m 



H , 



V,(z) 



/ 23^ 17 

+ (^12(1 + 6C2)C| + (-— + 12C2)C^Ca + Y^prif ) 

+ (-18C|, - 22C2 + 4C>/) In 



+ ( -^CICa - ^CfCI + Cluf + ^CpCAUf - \cpn) ] In 



+ 



162 



5 



CM 4 + I6C2 - 15C3 - 192C4 + — C2 + C'fCa -- - — C2 - 9OC3 + -C: 



17 326 



I^2(^ 

57 



6 9 

^^2^15503 188^ 1K2\ ^2 / 23 50^ 

+ ^.^.n,(-f^ + 6C2)+C.nj(g-lc.) 
+ (-6C3-Hc|C. + ^C^n,)ln3(||) 

+ {il - 32C2)C3 + - AC2)C',Ca - ^C^n,) In^ 
+ C^(6C. + fc.-in,)ln^(£|)ln(£|) 

+ CI{-9Cf - lie A + 2n/) In (^£f j In 
/ 21 

+ (^(- + 18(2-172(3)^3 + 



H , 



Mr 



,257 
108 



^(2 - K^)CICa + (-^ + f (2)C|n, ) In 



, m 



H y 



+ (^^^1 + "^-IcICa + ^^C,Cl - \n,Cl - ^-^CCau, + ^C^nj) In^ (^f" 



3 

, 289 2, ^ 10^ 9 



^4^2) CICa + - y(2)C^Cl + i(-5 + 8(2)C>, 



12 



54 



+ -C2)CFCAnf + ^C'i.n^ ) In 



m 



H , 



+ 



+ CfC\ 



297029 6139, 2509, 

+ ^(2 + ^(3 - 6(5 



1009 



C3 (-16(3 + 192(5 - 96(2(3) + CICa + ^(2 + ^(3 + ^(4 - 22(| - 23(2(3) 

187 



^Ci - y (2(3 



23328 324 ' 108 

421 47 179 40 19 2\ , . 

288 " 54^^ " l8^' ~ Y^' + y + ^^^^""^ 1 11664 " 



/31313 1837 



155 



23 
30' 



Cs + ^(2^ 



8 



+ 
+ 
+ 
+ 
+ 



+ 



729 



2V 27 
11 



121 



11 



21 



216 

.43 , C2 



54 

445 11 



m 



H , 



- - 9C2 + 38C3)C^ + (it + V - K^)CfCa + {-^- T7;C2)CpC 



108 12 



A 



7 (2 



289 C2 



+ i-^^ + ^i)CFCAnf + -C,n) In 



^4 



216 6 
121 



54 
11 



33 2^ 



In^l 
In 



41 

mij) 
( f^F 



In 



In^ 



143 53., 



■^ + yC2)C.Cl + (^-^C2)C>, 



+ ( 

-2 
+ 



289 C: 



108 3 

-8C2- 
, 245 



- -:^)CFCAnf - ^CfUj | In 



In 



m 



H , 



m 



H , 



33, 81^2N^3 .253 227^ 125^ 

yC3 + 96a--Cl)Ci + (- + -C. + ^C3 



57 

To 



(2)0 pC A 



67, 11 11 

\ C2 Cs 

48 18^ 12^ 10 



(2)^^(7^ + 



43 16 



144 9 



■C2 



14 



C3)C'Fnf 



,209 5, 7,,^^ 1^2\,/AiF\ 
+ (^-gC2 + -C3)C.C.n, + -C..j)ln^£|j 



+ 24(3^^ + 



7 



101 11 



415 



1^ 



+ y^^ + 12 



16 



C3 j C'pCa + 
178 22 



162 



9 - C2 - yC3 j C^Fnf + (— - yC2 - ^(3)^^^^^^^^, 



14 



+ (-- + -C2)C.n^)ln 



Vn(z) 



m 



(3.7) 



This is in agreement with Ref. 

The leading and some subleading colhnear corrections at NNNLO are 



a 



(3)C 
bb^H 

+ 
+ 



{-SCl ln^(l 



110 
~9~ 



CI^Ca 



20 



Cluf + 20C|, In 



^4 



m 



H , 



8CUl + 6C2) + C^fCa 



488 



+ 8C2 - 



ln^(l - z) 

^CfCI + Cf^{240Cf + UCa - 4n;) 



- 16^3 In^ 



+ Ci 



m 



H . 



16 



l2CF-^CA + ^nf]\n 



I^F 



m 



H , 



+ 



-l2Cl 



fclCA + Icluf ) In 



m 



H , 



ln3(l - 



bb->H at Tevatron S "=1.96 TeV ^1=111^^ 



gg->H at Tevatron S =1.96 TeV i^=m^ 



YOO 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 














NLOS/NLO 

NLOS+V/NLO 

NLO S+V+LC / NLO 

NLO S+V+NLC / NLO 




. . . 1 , , 


1 . . , 1 , , . 





160 



180 



(GeV) 



200 



= 0.8 



0.4 



0.2 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


















NLOS/NLO 






NLOS+V/NLO 






NLO S+V+LC /NLO 






NLO S+V+NLC / NLO 




, , , 1 . . , 


, , . 1 . , , 1 





140 160 
m„ (GeV) 



180 



200 



Figure 1: Left: The NLO ratios for bb —>■ H a.t the Tevatron. Here fi = fip = f^R = ^h- Right: 
The NLO ratios for gg H at the Tevatron. 
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Again we note that only the LC (ln^(l — z)) terms are complete. The NLC (ln'^(l — z)) and 
NNLC (ln^(l — z)) terms are not complete but, based on our study at NNLO, we expect them 
to be a very good approximation to the complete terms both in their analytical structure and 
numerically. Finally, in the ln^(l — z) terms we only show the scale terms that we expect to be 
exact at this accuracy. 



4 Cross sections for bb ^ H at the Tevatron and the LHC 

We now present a numerical study of the contribution of the corrections to the cross section for 
bb ^ H dX the Tevatron and the LHC. We use the bottom quark parton distribution functions 
(pdf) from the MRST2006 NNLO set of parton densities [23]. We are interested in the total 
cross section and the relative size of the higher-order contributions to it (the effect of parton 
radiation in transverse momentum distributions for this process has been studied in [23]). We 
also provide a few results for gg H using the gluon pdf from [23j to compare with bb —>■ H. 
In the results below we set the factorization and renormalization scales equal to each other and 
denote this common scale by fi. 

In Figure 1 we investigate the contribution of various terms to the complete NLO corrections 
for Higgs production at the Tevatron, with = 1.96 TeV and setting /x = mn, by plotting 
the corresponding ratios. It is important to clarify that in this figure NLO denotes the 0{as) 
corrections only (i.e. without the Born term). The left-hand side shows results for the process 
bb ^ H and the right-hand side shows for comparison results for gg ^ H. The curve marked 
NLO S / NLO denotes the percentage contribution of the NLO soft (S) corrections to the total 
NLO corrections. We see that this contribution does not surpass 50% for bb ^ H and 40% 
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Figure 2: Left: The NLO ratios for bb —>■ H at the LHC Here fi = fip = fJ'R = itt-h- Right: The 
NLO ratios for gg ^ H a.t the LHC. 



for gg H and thus the soft-gluon approximation is by itself inadequate. Adding the virtual 
terms to the soft, we get the soft plus virtual (S+V) approximation which, although better 
than the soft approximation alone, still does not provide a good approximation of the full 
corrections. Adding collinear corrections clearly substantially improves the situation. Simply 
by adding the leading collinear (LC) logarithms to the soft and virtual terms, the resulting 
S+V+LC approximation accounts for about 80% of the total NLO corrections for bb H. 
For gg H the S+V+LC approximation overestimates the total NLO corrections by a few 
percent. If we further add the next-to-leading collinear terms the approximation (S+V+NLC) 
gets even better for bb H, reaching over 85% of the total corrections. For gg ^ H the 
inclusion of NLC terms overestimates the cross section, but still by less than 10%. Clearly the 
inclusion of collinear terms greatly improves the approximation in both cases, and in particular 
for bb ^ H it is important to include the NLC terms. 

Figure 2 shows the corresponding NLO ratios for Higgs production at the LHC, ^/S = 14 
TeV. Again, the left-hand side shows results for the process bb H and the right-hand side 
shows for comparison results for gg H. We note that the results are quite similar to those for 
the Tevatron for both processes and our observations and conclusions are the same. In fact the 
soft and the S+V approximations are even more inadequate at the LHC (consistent with the 
fact that we are further away from threshold) which makes it even more essential to include the 
collinear corrections. We also note that in both Figure 1 and Figure 2 our results for gg ^ H 
are consistent with those published in [S]. Since gg ^ H has been studied at length elsewhere 
with results consistent with ours, we focus the rest of the presentation of our numerical results 
solely on the process bb — > H . 

In Figure 3 we investigate the contribution of various terms to the complete NNLO correc- 
tions for Higgs production via bb ^ H at the Tevatron (left-hand side) and the LHC (right-hand 
side), with /i = mn, by plotting the corresponding ratios. It is important to clarify that in 
this figure NNLO denotes the Ola"^) corrections only (i.e. without the Born term and NLO 
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Figure 3: Left: The NNLO ratios for bb —>■ H at the Tevatron. Here fi = fip = fJ'R 
Right: The NNLO ratios for bb ^ H at the LHC. 



niH. 



corrections). The curve marked NNLO S / NNLO denotes the percentage contribution of the 
NNLO soft corrections to the total NNLO corrections. We see that both at the Tevatron and 
the LHC the soft contribution is rather small, in fact even smaller than the relative contribution 
at NLO. The same holds for the S+V contribution. Inclusion of the leading coUinear logarithms 
improves the situation but is not by itself satisfactory since it only accounts for about 60% of 
the total NNLO corrections at both the Tevatron and the LHC. However, including the next-to- 
leading collinear logarithms vastly improves the approximation. The effect of the NLC terms is 
much more significant at NNLO than at NLO. As we saw in Fig. 1 and 2 at NLO the difference 
between the S+V+LC and the S+V+NLC curves was of the order of 5% at the Tevatron and 
10% at the LHC, so it was not overly significant. However at NNLO the the difference between 
the S+V+LC and the S+V+NLC curves is around 30% at the Tevatron and 40% at the LHC 
and thus the NLC terms are of utmost importance to get a good approximation. 

We actually plot two curves including next-to- leading collinear terms in Fig. 3. The curve 
S+V+NLC that we just discussed includes the full next-to-leading collinear logarithms. The 
curve S+V+NLCapp includes the approximate next-to-leading collinear logarithms from the 
expansion of the resummed cross section, Eq. (13. 6p . As noted before in Section 3, the expansion 
does not derive the full next-to-leading collinear logarithms; however, both at the analytical 
level and now as seen at the numerical level the difference between the two results is relatively 
small. At both the Tevatron and the LHC the S+V+NLCapp result accounts for around 
90% of the total NNLO corrections while the S+V+NLC result accounts for nearly 100% of 
them, which is rather remarkable. Finally, we also plot a curve (S+V+NNLC) that in addition 
includes the exact next-to-next-to-leading collinear terms (if instead we add the approximate 
NNLC terms from Eq. (13. 6p the resulting curve is practically indistinguishable because the 
NNLC approximate corrections are numerically very close to the exact NNLC corrections). 
From the figure we see that the NNLC terms alone do not make a large contribution, and that 
the S+V+NNLC results approximate the exact NNLO corrections very well. 
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Figure 4: Left: The K factors for bb ^ H at the Tevatron. Here fi = fip = f^R = ^h- Right: 
The K factors for bb H at the LHC 



Figure 4 shows the K factors, i.e. the ratios of the cross sections at various orders, for 
Higgs production via bb —>■ H at the Tevatron (left-hand side) and the LHC (right-hand side), 
with fi = rriH- Here and in the rest of the figures NLO cross section means the Born term 
plus the 0{as) corrections, NNLO cross section means the Born term plus the 0{as) and 
O(a^) corrections, and NNNLO cross section means the the Born term plus the 0{as) and 
C(a^) and 0{al) corrections. We note that we have used the same pdf for all curves because 
we are interested in the relative size of the terms in the perturbative expansion without the 
additional complications from different pdf. As the NLO / LO curve shows, the complete NLO 
corrections increase the LO result by around 60% at both the Tevatron and the LHC. The 
NNLO / LO curve shows that inclusion of the complete NNLO corrections futher increases the 
cross section by a substantial amount. The NNLO K factor is around 1.9 at the Tevatron and 
1.8 at the LHC. Finally, we include the approximate corrections (soft and coUinear) at NNNLO. 
We note that the soft corrections are complete and we plot one curve with the soft and leading 
collinear (S+LC) terms, another curve with the soft and approximate next-to-leading coUinear 
(S-|-NLCapp) terms, and a third with the soft and approximate next-to-next-to-leading collinear 
(S-|-NNLCapp) terms. We note that the difference between the S-I-LC and S-|-NLCapp curves 
is not very big, and between the S+NLCapp and S+NNLCapp curves it is quite small. Our 
investigation of the contributions of the soft and collinear terms at NLO and NNLO at both the 
Tevatron and the LHC gives us confidence that the NNNLO S+NNLCapp curves provide a good 
approximation of the complete NNNLO cross section. We note that the NNNLO S+NNLCapp 
K factor is between 2.00 and 2.08 at the Tevatron and between 1.86 and 1.97 at the LHC for 
Higgs masses ranging between 100 and 200 GeV. In both cases the NNNLO corrections provide 
a significant enhancement over the NNLO result. 

In Figure 5 we plot the cross sections for bb ^ H at the Tevatron (left-hand side) and 
the LHC (right-hand side). We show LO, NLO, NNLO, and NNNLO S+NNLCapp results 
for /i = TTiH. We note that the NNLO SV+NLC(app) or NNLO SV+NNLC(app) results 
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are indistinguishable from the exact NNLO curve on this plot. All five NNLO curves are 
on top of each other. This once again shows what we saw in more detail in figure 3, i.e. 
that our approximations including soft+virtual+next-to-(next-to)-leading coUinear terms are 
excellent. Also, as expected from Figure 4, the NNNLO S+NLCapp curve is indistinguishable 
from NNNLO S+NNLCapp. 

In Figure 6 we plot the cross sections for bb H at the Tevatron (left-hand side) and the 
LHC (right-hand side) and show NNLO and NNNLO S-|-NNLCapp results for two choices of 
scale, yU = mH/'2,2mH- We note that if we plotted instead the NNNLO S+NLCapp results 
they would be indistinguishable from the NNNLO S+NNLCapp curves. From the figure we see 
that the NNNLO S+NNLCapp result has very similar scale variation to the NNLO result. 

Finally, we calculate the uncertainty in the cross section from the parton distribution func- 
tions using the sets in [23] (pdf uncertainties for this process have also been studied in [25]). 
We find that the uncertainty is smaller at the LHC than at the Tevatron. For bb ^ H at the 
Tevatron we find that the pdf uncertainty varies from 3.1% for ttih = 100 GeV to 5.6% for 
tuh = 150 GeV to 8.0% for mfj = 200 GeV. For bb —>■ H a.t the LHC the pdf uncertainty varies 
from 2.0% for tuh = 100 GeV to 1.6% for tuh = 150 GeV to 1.3% for tuh = 200 GeV. We thus 
see that the pfd uncertainty is non-negligible and can be of the same order of magnitude as the 
scale uncertainty, especially for large Higgs masses at the Tevatron. 



5 Conclusions 

We have studied Higgs production at the Tevatron and the LHC via the channel bb H. 
We have calculated the complete soft corrections and the approximate next-to-next-to-leading 
collinear terms through NNNLO. We have shown that the inclusion of collinear corrections is 
essential in providing a good approximation to the complete cross section at NNLO. The soft 
and collinear NNNLO corrections provide significant enhancements to the cross section and 
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Figure 6: Left: The cross section for 66 —> if at the Tevatron with fi = = fiji = 2mH- 
Right: The cross section for hb ^ H dX the LHC 



must be taken into consideration for an improved theoretical prediction. The scale dependence 
of the cross section and the pdf uncertainties were also calculated. Analytical expressions for 
the soft and coUinear corrections through NNNLO were also provided for gg ^ H and the 
Drell-Yan process. 
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Appendix A: NNNLO soft and collinear corrections for 
Higgs and Drell-Yan processes: useful formulas 

In this Appendix we first provide detailed expressions for the quantities that appear in the 
resummed cross section, Eq. (12. ip . We then calculate a number of integrals for the soft and 
collinear corrections. Finally, we derive general formulas for the soft and collinear corrections 
for the processes under consideration in this paper. 



A.l Exponents in the resummed cross section 

For the quantity Ai in the first exponent of Eq. (12.11) we use the expansion Ai{as) = Af'^ 03/11 + 
Af'\as/7i)'^ + A^^\as/'n')^ + ■ ■ ■■ The results differ for quarks, i = q, and gluons, i = g. Here 
A« = Cf = {Nl - l)/(2iV,) with iVc = 3 the number of colors, and = Ca = N^- 
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= CfK/2 with K = (67/18 - C2) - 5n//9 [26], where Uf is the number of hght quark 
flavors {nj = 5 in our numerical results), and A^g^ = CaK/2; and Af^ is given by 
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2 /245 67 . 11 . 11 
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(A.l) 
(A.2) 
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Here and below (2 = 7rV6, Ca = 1.2020569 ■ ■ -, (4 = vr790, Cs = 1.0369278 ■ ■ -. 

Also Ui = (a,/7r)z/f ^ + (a,/7r)2z/f ^ + (a,/7r)3z/f ^ + ■ ■ -, with z/« = and z/« = Ca. It i 
convenient to lump together z/j with Ts as we discuss below. 

In the second exponent of Eq. (12.11) . ji/i is the moment-space anomalous dimension of the 
MS density [28l |29]. We write 



7^"^) In AT + 7, 
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and the parton anomalous dimensions 

7. = («./vr)7f ^ + {as/i^flf^ + ■■■ 
with 7« = 3C^/4, 7« = /5o/4, 



and 
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The (3 function in the third exponent of Eq. (12. ip is given by 
1 da^ 



lag dlnfi 
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where g'^ = Arras, with /?o = (HCa - 2n/)/3 and /?i = 34^2/3 - 2n/(C^ + 5Ca/3) []. Note 
that 
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Also da^ is for bb ^ H and qq V and 2 for gg H. 

The hard and soft functions can be expanded as if = a^"'" H^^^+{a'^"=^^ /7r)H'^^^ + {a'^"''~^'^ /tt'^)!!'^'^'' 
+ (af-+V7r3)ii(3) + ■ ■ ■ and ^ = 5^ + (a,/7r)5(i) + {ajTrfS^^^ + (a./Tr)^^^^) + .... At lowest 
order, the Born cross section for the partonic process is = a^"" H^'^^ S^^\ 

The soft anomalous dimension is denoted by and can be expanded as r5 = (as/7r)r^^^ 
+ {as/'n')^Tg'^ + {as/TT)'^T[^^ + ■ ■ ■■ For the processes bb ^ H and qq V, ^^sqq — ^f- For the 
process gg H, F^^^^ = Ca- If we define G^^^ = F^"-* — i/f"^ we can get explicit expressions 
through [TSl EOl EI]. We have C'gj = = 0, 
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Gj? = ^G^S ■ (A.13) 



A. 2 Useful integrals 



In the resummed cross section we encounter certain integrals that need to be evaluated. They 
are of the form 





"ln"(l-2)' 
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The expressions for J„ have been presented up to n = 7 in [21] and including 1/A^ terms up 
to n = 3 in [12\. Here we extend these results by including 1/A^ terms up to n = 5, which are 
needed for NNNLO expansions: 
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h{N) = In^ iV - 2C2 In^' iV - 4C3 ln2iV-3(C2' + 2C4) In iV- 4(^(2(3 + ^(5 
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We also encounter integrals of the form 

J^{N) = dzz^-Hn''{l-z). (A.16) 
Jo 

The expressions for J„ have been presented up to tt, = 3 in [12]. Here we extend these results 
up to n = 5, which are needed for NNNLO expansions: 
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A. 3 General formulas for soft and collinear corrections 



Below, we write general expressions for the soft and collinear corrections at NLO, NNLO, and 
NNNLO for the processes under study. Detailed results are provided in Section 3 for bb H, 
in Appendix B for the Drell-Yan process, and in Appendix C for gg H. We note that for 
processes involving a quark running mass, such as bb H, the expressions below are valid for 
the running mass at scale M. If the mass is evaluated at scale fiji (as in Section 3) we have to 
include additional terms as discussed at the end of this appendix. 
The NLO soft and virtual corrections are given by 
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with z — js. If we define C, to denote Cp for quarks and Ca for gluons then the coefficients 
are given by 

C3 = 4C, (A.20) 
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(A.23) 



and Ti the scale-independent virtual corrections for each process. We note that only involves 
the scale dependence. This simplifies the structure of the results at higher orders. 
The NLO coUinear corrections are given by 
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where c\ = —ACiKi and = 2CjKj ln(/i|,/M^) + 2Ci, where Kg = 1 and = 2. 
The NNLO soft and virtual corrections are given by 
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where V^'^^ denotes the two-loop virtual scale-independent terms for each process, which cannot 
be derived from resummation. 

The NNLO (approximate) coUinear corrections are given by 
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where T2 = 2Cj denotes the scale-independent part of C2. Note that only the leading terms are 
complete. However, the ln^(l — z) and ln(l — z) are nearly complete as discussed elsewhere in 
the paper. We do not calculate the constants other than the scale-dependent terms as described 
elsewhere in the paper. 

The NNNLO soft-gluon corrections are given by 
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The NNNLO (approximate) coUinear corrections are given by 
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Note that only the leading terms are complete. However, the ln^(l — z) and ln'^(l — z) are 
expected to be nearly complete as discussed elsewhere in the paper. We do not calculate 
ln^(l — z) or lower terms other than some scale-dependent terms as described elsewhere in the 
paper. 

In the calculation of 66 if we also have to consider the renormalization scale logarithms 
from the MS bottom quark running mass, m^. Since n'^drrib/diJ? = —•jm^b with 7^ the mass 
anomalous dimension 
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Appendix B: Soft and collinear corrections for the Drell- 
Yan process through NNNLO 

The complete NLO soft and virtual corrections are 
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where F^y is the Born term, whose exact form depends on which gauge boson is involved in 
the process. 

The leading and next-to-leading collinear corrections at NLO are 
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The complete NNLO soft and virtual corrections are 
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This in agreement with Refs. [7l [33] where the virtual terms were calculated in [33] 
The leading and some subleading coUinear corrections at NNLO are 
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Several remarks are in order here. The leading coUinear (LC) logarithms (i.e. ln^(l — z)) are 
complete. The next-to-leading coUinear (NLC) logarithms (i.e ln^(l — z)) are not complete. 
However, numerically they are an excellent approximation to the complete NLC terms (84% at 
fiF = Q)- Also analytically the CfCa, ujCf and the ln(/Xp/(5^) terms are exact. The difference 
of exact minus approximate NLC terms is 7C|^/2. The next-to- next-to-leading coUinear (NNLC) 
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logarithms (i.e ln(l — z)) are also not complete. However, again numerically they are an 
excellent approximation to the complete NNLC terms (107% at fip = Q). Also analytically 
the CfCaC-i. In^ifil/Q^): and ln(/i|/g2) terms are exact. The difference of exact minus 

approximate NNLC terms is 7C|,/4 - poCp + 5CfCa/4 - 2C'^ln{nl/Q^). Finally, in the 
constant terms we only show the scale-dependent terms; the In^ {/ip / Q'^) , ln(/i|,/(5^) ln(/i|j/(5^) 
and ln(/i|j/(5^) terms are exact while the ln(//|./(5^) terms are almost exact. We note that we 
can even get some additional constant terms such as —IGCpCs- 
The complete NNNLO soft-gluon corrections are 
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The scale-independent terms in the above expression are in agreement with Refs. [T8J [19] 
The leading and some subleading coUinear corrections at NNNLO are 
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Again wc note that only the LC terms (i.e. ln^(l — z)) are complete. The NLC (i.e. ln^(l — z)) 
and NNLC (i.e. ln^(l — z)) terms are not complete but, based on our study at NNLO, we 
expect them to be a very good approximation to the complete terms. In the ln^(l — z) terms 
we only show the scale terms that we expect to be exact at this accuracy. 



Appendix C: Soft and collinear corrections for gg 
through NNNLO 



The NLO soft and virtual corrections are 
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where v = 246 GeV is the Higgs boson vacuum expectation value. 
The leading and next-to-leading collinear corrections at NLO are 
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The complete NNLO soft and virtual corrections are 
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where is the top quark mass. This is in agreement with [Sj HI [5], [3l], and the virtual 
corrections were calculated in [35] . 

The leading and some subleading coUinear corrections at NNLO are 
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Again, only the leading coUinear (LC) logarithms (i.e. ln^(l — z)) are complete. The next- 
to-leading collinear (NLC) logarithms (i.e ln^(l — z)) are not complete. However, numerically 
they are an excellent approximation to the complete NLC terms (83% at fip = itlh). Also 
analytically the ujCa and the ln(/i|,/m|^) terms are exact. The difference of exact minus 
approximate NLC terms is 7C\/2. 

The NNNLO soft-gluon corrections are 
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The scale-independent terms in the above expression are in agreement with Refs. 
The leading and some subleading coUinear corrections at NNNLO are 
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Again, only the leading terms in the above expression are exact. 
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